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Abstract 

The paper compares computational aspects of four approaches to compute conser- 
vation laws of single differential equations (DEs) or systems of them, ODEs and PDEs. 
The only restriction, required by two of the four corresponding computer algebra pro- 
grams, is that each DE has to be solvable for a leading derivative. Extra constraints for 
the conservation laws can be specified. Examples include new conservation laws that 
are non-polynomial in the functions, that have an explicit variable dependence and 
families of conservation laws involving arbitrary functions. The following equations are 
investigated in examples: Ito, Liouville, Burgers, Kadomtsev-Petviashvili, Karney-Sen- 
Chu-Verheest, Boussinesq, Tzetzeica, Benney. 

1 Introduction 

As is well known, conservation laws play an important role in mathematical physics. The 
knowledge of conservation laws is useful in the numerical integration of partial differential 
equations (PDEs) [20], for example, to control numerical errors. Also, the investigation of 
conservation laws of the Korteweg de Vries equation was the starting point of the discovery of 
a number of techniques to solve evolutionary equations [2H] (Miura transformation, Lax pair, 
inverse scattering technique, bi-Hamiltonian structures). The existence of a large number of 
conservation laws of a PDE (system) is a strong indication of its integrability. Conservation 
laws play an important role in the theory of non-classical transformations [221, and in 
the theory of normal forms and asymptotic integrability |24j . Programs described below are 
able to find conservation laws involving the independent variables explicitly. Finding such 
conservation laws is a good challenge for the inverse scattering technique. 

The purpose of the methods described below is to pose as few restrictions as possible 
on the differential equations (DEs) to be investigated. For example, it is not assumed that 
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any Lie symmetries are known, nor that the equations are equivalent to the Euler-Lagrange 
equations of a variational problem. Instead we attempt to solve the conservation law condition 
directly. The strategy will be to make a local ansatz involving only the dependent variables 
and their derivatives. Further, the order of the derivatives is bounded in order to obtain 
an over-determined PDE problem which subsequently is solved with the computer algebra 
package Crack [33]. |33|. 

In an earlier paper (35] three of the methods were discussed with emphasis put on the 
computer algebra algorithms involved. In this paper we present an additional fourth method 
and compare these methods in terms of complexity and functionality. 

The rest of the paper is organised as follows. After specifying the notation that is used, 
in section |HJ a reminder on issues of the equivalence of conservation laws will provide the 
motivation for the four approaches which are explained in section 0] followed by an overview. 
In section |S] related computer algebra programs are shortly described and examples are given. 
Extensions of the basic usage of these programs are discussed in section 

2 Notation and setup 

We adopt the notation of the book of Olver |27| where the question of equivalence of con- 
servation laws is described in more detail in Chapter 4.3. This compact notation will be 
explained using the sine-Gordon equation which will be used as an example throughout the 
paper. 

• In general derivations or generally applicable formulas the independent variables are 

p ). In examples independent variables are t,x and y. 

• Functions are denoted in general by u = (u 1 , u 2 , . . . , u q ). In examples the function is u 
and in a few examples an additional function is v. 

• Partial derivatives are written as lower index as in u tx — sin u — 0. If partial derivatives 
are repeated n times (n > 2) then this may be indicated by writing n in front of 
the variable in the index like in d 7 u/(dt 4 dx 3 ) = u^x (which is especially used in the 
appendix due to the high derivatives occuring there). 

• An upper index in brackets like denotes the set of all derivatives of all components 
u 1 , . . . ,u q of order up to n including order zero, i.e. the functions themselves. In our 
example with u = u(t,x) this would, for example, mean = {u, u t , u x , u tt , u tx , u xx ] . 

• The differential equations that are to be investigated concerning conservation laws are 
denoted = A(x, u^), in our example A(x,u^) = u tx — sinu. If we have q functions 



2 



u 1 , . . . , u q then it is assumed that a system of q equations = A 1; = A 2 , . . . , = A q 
is given. 

• Latin indices i,j have a range 1, . . . ,p used for x % and Greek indices a, \i have a range 
1, . . . , q used for u a and A M . 

• Whereas u^> denotes the set of all possible indices up to order n we will also need a 
way to specify any particular partial derivative of an unspecified order. For that we 
will use the so-called multiple index j. Then any u x ,u t ,U2t5 X are all examples for uj. 
An example for the use of j is the notation of total derivatives: 

• For the total derivative we use the symbol D: 



• If a relation = W(x, «(")) is said to be satisfied identically for any solutions of the 
equation(s) = A then this means that W can be written as a linear combination of 
A and any total derivatives of A with arbitrary x— and u^ — dependent coefficients 
which are non-singular for = A. The notation adopted for that relation is 



which may be described as U W vanishes modulo A = 0". The way to check this 
circumstance in practice is 

— to solve each A M for one leading derivative (such that no other derivative in A M is 
a derivative of its leading derivative and that no two leading derivatives from two 
different A M , A u coincide and that none is a derivative of another one), 

— to substitute in W(x, u^) all these eliminated leading derivatives as well as all 
derivatives of the leading derivatives. 

— We have = W| A=0 if and only if W vanishes identically after all possible substi- 
tutions have been performed. 




Using the multi index j a more compact notation is: 



D x i = d x i + Y,Y, u Jx^ a j- 



a=l J 
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Example: For the sine-Gordon equation = A = u tx — sin u we will take u tx as leading 
derivative and do substitutions u tx = sinu. Let us assume an expression 

W = —8u xx u xxt + 4:U x u tx + 8u x u xx cos u — 4u x sin u (1) 

is given and we want to compute W^|a=o- The highest derivative of u tx is u xxt = d x u tx . 
We therefore substitute u txx = d x (smu) = u x cosu giving 

W = 4u 3 x u tx - Au x sin u. (2) 

Substituting now u tx = sinu gives W = 0. We therefore found that W as given in ((TJ) 
satisfies M^|a=o = 0. 

A conservation law of A will be given in form of a so-called conserved current P l where 
each P % is a differential expression in x and u, i.e. P % = P l (x, u^) which has a vanishing 
divergence due to = A. In other words 



0= Y, D * iPl 



= DivP| 

A=0 



A=0 



Example: For the current 

P» = -4u^ + ttl, P x = 4u 2 xC osu, (3) 

we find 

DivP = D t F* + D X P X = -8u xx u txx + Au 3 x u tx + 8u x u xx cos u - Au 3 x sin u, (4) 

which is equal W in ((TJ). Because of D1vP|a=o = W^|a=o = the vector P in (j3J) is a 
conserved current of the sine-Gordon equation and represents a conservation law. 

About the conserved quantity: By integrating = DivP over a region of the p- 
dimensional x l -space we obtain the vanishing of a surface integral over that region: 

= fpidSt. 



By choosing a region with a cylinder-like shape where the radius of the cylinder is very 
large and the bottom and top side of the cylinder lie in the p — 1 dimensional surfaces 
x 1 = a = constant, and x 1 = b = constant the surface integral takes the form 

= / P l dx 2 . . . dx p - I P 1 dx 2 ...dx p 



if we assume that P 1 falls off sufficiently quickly for any x l — > oo and therefore the 
integral over the mantle of the cylinder vanishes S maxA \ c P l dSi — > 0. The minus sign 
comes in because the normal vectors of the top and bottom surface point into opposite 
directions. 

If one of the coordinates plays the role of time, say x 1 , then / P 1 dx 2 . . . dx p is called a 
constant of motion as it does not change from any time a to any time b. 

Example: Continuing our example we have as a constant of motion 

/oo 
{-Au x l + ul) dx 
-oo 

if u is a solution of the sine-Gordon equation. 



3 The equivalence of conservation laws 

Although two conservation laws may look rather different, i.e. their two conserved currents, 
say P and P, may be different, nevertheless their information content may be the same. To 
have a method of counting conservation laws and of comparing them, we need a unique way 
to characterise them. To do that we first look at ways conservation laws can look different 
but be equivalent. Two conservation laws DivP = and DivP = are equivalent if 
= Div(P — P) — DivP is a trivial conservation law. 

(i) The first kind of equivalence of two conservation laws is the case that R = for all 
solutions of A = 0, i.e. if P and P differ only by multiples of A and by total derivatives of 
A (i.e. by DjA). To test whether this is the case for two given conserved currents P and P 
one has to check whether 

= (P - P) 

v ' A=0 

holds. As described above one solves the equation = A (or system of equations = A M ) for 
the leading derivative (s) and substitutes that (them) in P — P. If P — P becomes identically 
zero then the two conservation laws based on P and P are equivalent. 

If the conservation laws are not yet calculated and one wants to ensure that the compu- 
tation of P gives a unique result, without arbitrariness due to terms vanishing because of 
A = then there is no need to solve A = for some leading derivative(s) uj. In that case 
one just drops from the beginning of the calculation the dependency of P on the leading 
derivative(s) uj and all derivatives of u j. 

Example: When computing conserved currents P for the sine-Gordon equation = u tx — 
smu, then in the ansatz for P the components P*, P x are assumed to be independent of u tx 
and derivatives of u tx , like u ttx ,u txx , ■ ■ ■ (see the appendix). 
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(ii) The second kind of equivalence of two conservation laws occurs if R l = P % — P l = 
J2j DjV' 1 ^ for some expressions V i ^(x,u (jn ') = —V^ 1 , anti-symmetric in i,j, because in that 
case DivP = YU,j DiDjV^ = (due to the symmetry DiDj = DjDi and the anti-symmetry 
y % i = —V^) for any functions u, not necessarily solutions u(x) of A = 0. The existence of 
V ij (x, u^) satisfying P i = P*+Ej DjV ij may not be obvious and may require a computation 
checking DivP = DivP. For ODEs this problems does not occur as there is only one 
independent variable and no antisymmetric . 

The solution to this problem is not to compare conservation laws by comparing their 
conserved currents P and P but by comparing them by their integrating factors, for PDEs 
they are called characteristic functions, in the following way. For a conservation law to 
satisfy DivP| A=0 = means that DivP is identical to a linear combination of A^ and 
total derivatives DjA^. Partial integration can rewrite that as a divergence plus a linear 
combination of the A M alone: 

DivP| A=0 = (5) 
<^=> 3Qi : DivP = Q J U D jA u (identically in all x, u a j) (6) 

U,J 

= ^2 Dj(Q^A„) — Dj(Q J u )A u (repeated partial integration) 

U,J 

= DivP + ]T Q v A v 

V 

DivP = Y.Q" A » (after renaming (P-P) ^ P). (7) 

The integrating factors Q u are called characteristic functions as it is known (|27j. p. 272) 
that for a totally non degenerate system A = 0, the equivalence class of conservation laws 
DivP|A=o = is characterised uniquely by the functions Q u up to equivalence of type (i). 

One can look at equation (J7J) as a determining equation for P and Q v as in the method 
described in section below. Alternatively one can formulate a system of conditions that is 
equivalent to (JJJ) but which involves only functions Q v . That is achieved using the property 
of Euler operators (also called variational derivatives) E u = J2j{~D)jd/duj when acting on 
an expression they give identically zero iff this expression is a divergence. Conditions for the 
Q v are therefore 

V^: = E v (V Q,A,) = Y,{-D)j V Q,a\ . (8) 



Example: Allowing Q to be of 2 nd order, i.e. Q = Q(t,x,u,u t ,u x ,Uu,u xx ) where we dropped 
the dependence on ut x which is equal to sin u (see remark above) we find that 

Q ■ A — Q • [u tx - sin u) 
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does depend on u } u t} u X} u t t,u tx ,u xx and the Euler operator therefore reads 

E = av-D^-Dtd^ + i-DtX-Djd^ + i-DJi-Djd^ + i-DJi-Djd^ 
= d u - D t d ut - D x d Ux + D 2 t d utt + D t D x d Utx + D% xx . 

Requiring condition(s) (JHJ) to be satisfied identically in all x\ u a and derivatives of u a (i.e. 
uj) is equivalent to the condition (JJJ). System (jHJ) is often very large. It can be considerably 
shortened if it is projected onto the space of solutions |a=o ( as described above): 

Vz/. (9) 

A=0 

Conditions @ are known as adjoint symmetry conditions which are necessary but not suffi- 
cient for the Q M to be characteristic functions of first integrals. 

iii) For any two conservation laws = Div P and = Div P, = Div(P + P) is also a 
conservation law. By determining conservation laws with characteristic functions of succes- 
sively increasing order, constant multiples of characteristic functions of lower order can be 
dropped. 

iv) In the case of (systems of) ODEs the characteristic functions are called integrating 
factors, and P is a scalar, called a first integral. Any arbitrary function of first integrals is a 
first integral as well. 

The four approaches described in the following four sections are to solve conditions (0), 
0, © and p. 

4 The four approaches 

4.1 A first approach 

The first approach is to solve 

BivP\ A=0 = (10) 

directly. 

The condition (fTUj) is made over-determined by restricting the P l to be differential ex- 
pressions in the u of at most some order k, i.e. P 1 = P l (x,u^ k '). Characteristic features of 
this approach are 

(+) A single, first order PDE involving only few terms is to be solved. 
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(0) Characteristic functions have to be computed from P in a straightforward calculation 
(described in jUS])- This is done within the computer algebra program ConLawI 
which implements the first approach including the computation of related characteristic 
functions Q^. 

(— ) It would be computationally expensive for a corresponding computer program to drop 
during the process of solving (jl(Jj) any free functions V 1 ^ = —V^ 1 (see the discussion 
of the second kind of equivalence of conservation laws in the previous section) which 
correspond to trivial conservation laws 1 . Hence, the condition (fTUj) has to be solved 
first in full generality and trivial conservation laws (i.e. V lJ ) have to be identified and 
dropped afterwards. That means that the task for the computer program is made 
unnecessarily hard by the presence of the trivial conservation laws in the general solution 
of (jl(J|) . A rule of thumb says that the difficulty in solving a linear over determined 
PDE system depends less on the order or size of the PDE but more on the complexity 
of the result 2 . That means the trivial conservation laws will complicate the solution of 
()1U|). the more so the more independent variables are present. 

(— ) In most cases the expressions for the P l are more complicated than the expressions 
for the characteristic functions which by the above rule of thumb indicates a more 
difficult computation than the solution of equations involving only Q^. 

To illustrate and compare all four approaches we will apply each to finding conservation laws 
of the sine- Gordon equation 

u tx - sin(u) = 0. (11) 

If the program ConLawI is called to find conservation laws with conserved current P*, P x 
of order 0, then it will reply that it is not applicable. This is because Div P would be of first 
order in u, so equation (|TT|) could not be used to substitute u tx (when computing |a=o in 
(fTU|)) and therefore any conservation laws found would be valid for any function u(t,x), not 
necessarily only for solutions of the sine-Gordon equation. These conservation laws would 
therefore be trivial, falling into category (ii) in section El 

Details of higher order investigations are given in table 1 below, u^ n > stands for all 
derivatives of u of order to n. u tx n ^ stands for all derivatives of u tx up to order n, for example, 
Ufa would be the derivatives u tx ,Utt x ,u txx . Finally, u^/u^ stands for all derivatives of u 

1 An algorithm for that is given in |35| . 

2 For example, if an over-determined PDE (system) has no solution then a differential Grobner Basis 
calculation will quickly produce PDEs of lower and lower order until a contradiction is reached. On the other 
hand, if a PDE system has arbitrary functions in its general solution (as is the case with the PDE l|10[l ) then 
computing a differential Grobner Basis will not produce a system that is solvable by only integrating ODEs, 
it will involve PDEs. 
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up to order n apart from u tx and all its derivatives up to order k. The conservation laws are 
given in the appendix (in the table only the equation number is cited). For each conservation 
law in the appendix (apart from the first) there exists another one resulting from the exchange 
t x. 

The times given in the table are measured on a 266 MHz Pentium PC running a 80 MByte 
Reduce 3.6 session under Linux using the Sep. 1998 version of the program Crack for 
solving the over-determined conditions. The 80 MByte were not necessary. For example, it is 
possible (using ConLaw2 which implements the 4 th method described below) to investigate 
up to 4 th order laws with 4 MByte and up to 7 th order laws with 8 MByte. To get this high 
in order with relatively low memory consumption, one has to give in Crack the study of 
integrability conditions a higher priority than the integration of equations. The price is a 
higher computing time. The times in the last column are to be understood only as very rough 
indicators 3 . They depend sensitively on the order of priorities with which modules are to be 
used within the program Crack (see the manual [33] and about its availability the end of 
the section |7j). 

When condition (jlOj) is solved, the P l that are computed initially do not contain u tx nor 
its derivatives. Afterwards a computation as outlined in - (J7J) is performed such that 
finally ConLawI is able to return the conservation law in the form (J7J). In the process of 
computing this form (JJJ) the new P l may now involve Ut x (through R in © - ©)• 



order 
of P i 


no of 
terms 


independent 
variables, [no of var.] 


functions to 
compute, [no of arg.] 


cons, 
laws 
found 


time to 

solve jini) 


1 


8 


t } x }U W/u tXl [7] 


P\P x (t,x,u { V), [5] 




9 sec 


2 


12 


t,x, U W/u£\ [9] 


P\P x {t,x,u^/u tx ), [7] 


(ED 


38 sec 


3 


18 


t,x,vF>/ug\ [11] 


P^P^x^^/u^), [9] 


none 4 




4 


26 


t,x,u^/u tx z \ [13] 


P\P*{t,xM 4) /u tx 2) ), [11] 


none 5 





3 For example, the computing times reported in are at the time of revision of this paper (July 1999) 
already reduced by a factor of more than ten for higher orders. 
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Table 1: The program ConLawI applied to compute conservation laws of the 
sine-Gordon equation. 



4.2 A second approach 

The next approach consists in solving 

Div P = ^Q V A V (12) 

V 

directly, i.e. finding P*,Q M that satisfy (fT2"j) identically in x l ,Uj. Equations A = DjA = 
are not used for substitutions in (fT2J) but they are used to reduce dependencies of the Q^. 

The problem becomes over-determined by restricting the order of the Q^, i.e. = 
Qn(x,vS k >) for some k and by taking <— C^|a=o, i-e. having independent of one 
leading w-derivative (and their derivatives) from each one of the equations A v . If would 
be allowed to depend on all u^ 1 ' which occur in (JT2J) then this equation could simply be solved 
algebraically, by eliminating one of the Q^. But that would mean division through one A M 
and therefore being singular for solutions of A^ = 0. 

The second approach has the following characteristics: 

(+) The conservation law condition (|12|) is a single first order PDE as in the first approach. 

(+) By calculating characteristic functions and furthermore characteristic functions 
Qh\a=o, conservation laws are uniquely characterized. 

(+) The effort in formulating conditions is as low as in the first approach. 

(0) The P % and are computed in one computation. 

(0) The number of functions to compute is higher than in the first approach and also the 
number of derivatives of u on which these functions depend on because no substitutions 
are done in (|12p. The resulting complication is not too big as more variables means a 
higher over determination and simplification. 

4 Crack was not able to solve all the equations completely because the general solution of I jlOjl involves free 
functions (related to trivial conservation laws) which complicates the problem considerably for the computer 
program. 

5 The computer memory was not sufficient to complete the computation. 



10 



If the order of A is n and the order of is chosen to be k then the order of P % at 
the start of the computation can be assumed without loss of generality to be max(A;, n) 
(see |26j). 'Without loss of generality' means that a trivially conserved current P can 
be subtracted from P such that P — P is of order max(fc, n). If the right-hand side of 
equation (fT2j) is known to be linear in the highest derivatives of order max(fc, n) then 
P % at the start of the computation can even be assumed without loss of generality to 
be of order max(fc, n) — 1. 

In this approach the investigations with k < n are not much simpler than the case 
k = n. This matters when the order n of A and the number p of variables x are high. 
Therefore this approach is not very efficient for low order conservation laws of high 
order equations. 

For example, for zeroth order conservation laws (k = 0) of the Kadomtsev-Petviashvili 
equation (fTTJj) the P l are taken initially as functions of the 23 variables t, x, y, u, u t , 
u x ,u y , u u , ■ ■ ■ , u yy , u m , • • • , Uyyy and the conservation law condition (fT2J) is a condition 
in 38 variables (including the 4 th order w-derivatives). That is a much harder problem 
than the corresponding conditions (JEJ) , © ■ For example, in this case condition (JSJ) is 
a single 4 th order PDE in also 38 variables but for only one function Q of only four 
variables! 

When looking for conservation laws with the first method, gradually increasing the 
order of the conserved current P gives each conservation law in its lowest order form, 
i.e. a form where P is of minimal order. This is not necessarily the case using the 
2 nd method. The transformation ©-0 adding R to P may increase the order of P. 
This implies an increase of complexity having to go up in order to get the equivalent 
conservation law. To give an example, the Tzetzeica equation u x t = e u — e~ 2u (analysed 
in [HO], [21]) nas the conservation law 



= D t 



*^^XXX ^^xx ^^^xx^x ^x 



+ 

-3e« (u 2 xx + 2u xx u 2 x + 2ut) - 3e~ 2u (2u 2 x - 8u xx u 2 x + - 1 



u,. 



with a third order conserved current. (In [21] an infinite list of conservation laws is 
given.) Bringing the above conservation law to the form (fT2"j) as it would be found with 
the second method, it becomes 



6 ^Uxxxxx ~\~ ^^xxx^xx Stt^xx^x ^ii"^ ^x) i^tx & ~\~ & 



■2u 



D 



3^xxx ^^xx ^-^^xx^x ^x 



+ 



D x 3 2,Ut x u xxxx 2,Ut xx ii xxx -\- 5ut x ii xx I0ut x u xx u x 
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-\-e \2,u xxx u x ^u xxxx 6it xa , -|- 8u xx ii x 1u x j 

-\~€ ( 2u&xxx ^^xxx^x 3u xx 2,U XX U X U x 

with a 4 th order conserved current. 



Applying the program ConLaw3 that corresponds to the above method to the sine-Gordon 
equation (JTTJl gives the following table. 



order 
of Q 


no of 
terms 


independent 
variables, [no of var.] 


functions to 
compute, [no of arg.] 


cons, 
laws 
found 


time to 
solve jnj) 





10 


t,x,u( 2 \ [8] 


P\P x (t,x,u^), [5] 
Q(t,x,u), [3] 


none 


3 sec 


1 


10 


t,x,u( 2 \ [8] 


P\P x ,Q(t,x,uW), [5] 




8.3 sec 


2 


10 


t,x,u( 2 \ [8] 


P t ,P x (t,x,u^), [8] 
Q(t,x lU W/u tx ), [7] 


none 


3.5 sec 


3 


16 


t,x,u&\ [12] 


P t ,P x (t,x,u ( -^), [12] 
Q^x^u^/u^), [9] 


none 6 





Table 2: The program ConLaw3 applied to compute conservation laws of the 
sine-Gordon equation. 



4.3 A third approach 

Instead of calculating the conserved current P % directly, the third approach is to calculate 
characteristic functions first (see, e.g. Proposition 5.33 in [27]) and from them P l af- 
terwards using formulas of Anco & Bluman 0,0,11] i n a form described in j^J] or using 
repeatedly the Crack routine for integrating exact DEs (see also the section on Homotopy 
Operators in (23 )• The condition (JHJ) (as derived in [2Zj,|2]) is: 

= E(-^(t^£^ A m) V ^ (13) 
6 The computer memory was not sufficient to complete the computation. 
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Typical features are: 

(+) Equations (fT3"j) are equivalent to (fT2*|l and therefore necessary and sufficient. 

(+) The usually more complicated P l are eliminated and as in the 2 nd method, no trivial 
conservation laws are calculated which otherwise unnecessarily complicate the calcula- 
tion. 

(+) The highest w-derivatives in conditions ()13|) are of the order 2n where n is the order 
of the w- derivatives in ^Q^. The harder the problem, i.e. the higher n and the 
higher the number of variables, the more w-derivatives occur only explicitly in (|T3"j) and 
can be used for a direct separation (splitting). Higher over determination simplifies the 
solution of (|T3|). 

(— ) Equations (JT3J) consist of as many equations as there are dependent variables u^. The 
unknown functions appear with n th order derivatives. 

(— ) For an increasing order of the Q^, increasing number of u v and increasing number of 
x l , the size of (II 3j) can soon become unmanageable. 

Applying the program ConLaw4 that corresponds to the above method to the sine-Gordon 
equation gives the following table 3. The striking feature of this approach is the quick 
increase of the "size of conditions" . Apart from the order case they increase by a factor 
of about 7 which itself is increasing slightly with the order. The size of conditions prevents 
going higher in the order. On the other hand, the completeness of the generated conditions 
simplifies the solution in difficult cases and speeds up the solution of the over determined 
system as long as it is not already too large at the beginning. 
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order 


no oi 
terms 


independent 
variables, [no of var.] 


functions to 
compute, [no of arg.] 


cons, 
laws 
found 


time to 
solve (HHD 
h:min:sec 





7 


t,X,U^,U tx , [6] 


Q(t,x,u), [3] 


none 


0.7 sec 


1 


22 


t,x,u (2 \ [8] 


Q(t,x,uW), [5] 




2.8 sec 


2 


154 


t,x,u®, [17] 


g(t,x, M ( 2 )/^), [7] 


none 


4.7 sec 


3 


1116 


t,x,u^\ [24] 


Q{t,x,uW/u£>), [9] 


iU 


5 min 17 sec 


4 


8402 


t,x,u^ 5 \ [34] 


Q(t,x, M ( 4 )/ M i 2) ), [11] 


none 


10 h 49 min 7 


5 


64064 


t,x,u (6 \ [41] 


Q(t,x,u^/u t ^), [13] 




> 2 days 



Table 3: The program ConLaw4 applied to compute conservation laws of the 
sine-Gordon equation. 



4.4 A fourth approach 

Projecting conditions (fT3j) into the space of solutions of A = we obtain 



Vis. 



(14) 



A=0 



The characteristic features of this method are similar to those of the third method with the 
following modifications: 



7 This time was nearly completely spent to formulate the condition and to separate it into 823 individual 
equations for Q. Then already the 3 rd step gave that Q cannot depend on 4 th order derivatives. 
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(+) The conditions usually involve fewer terms than in the third approach which can be 
decisive. But as the conditions (|14|) are not sufficient, they are less over determined 
and may be harder to solve than those in the third approach. 

(— ) Because the conditions (JHJ) investigated by this fourth method are not equivalent to 
and sometimes less restrictive than conditions (fT3|) . the solutions of (fT4*Jl need not 
represent conservation laws. Therefore after computing the from (jl4j) . it has to be 
checked whether P i exist that satisfy DivP = X^QjA, (0, 0,11, El)- If they do not 
exist then the Q M correspond to an adjoined symmetry but not to a conservation law. 

(— ) If the fourth method finds adjoined symmetries which are not conservation laws, then 
it is still possible that these adjoined symmetries can be combined linearly to give 
conservation laws. But how to combine adjoined symmetries to give conservation laws 
is not answered by solving (|14|) . it has to be investigated separately. This will be 
illustrated with the following simple example. 

Applying the third method through the program ConLaw4 to the ODE u" + u = 
and restricting the search to integrating factors Q = a(x)u + b(x) that are linear in u, 
the program finds the following 5 integrating factors 

cos(x) 2 u' + cos(x) sm(x)u, 2 cos(x) 2 u — 2 cos(x)u' sin(a;) — u, cos(x), sin(x), u', 

and therefore it finds 5 first integrals. In comparison, using the same ansatz for Q but 
now using the fourth method with the less restrictive condition (fl"4"j) the corresponding 
program ConLaw2 finds 8 solutions for Q: 

—2cos(x) 2 u' — 2cos(x)sm(x)u + u', cos(x), sin(a;), u', 

cos(x)u'u + sin(x)-u 2 , cos(x)u 2 — u' sin(a;)M, — cos(x) 2 u + cos(x)v! sin(a;), u. 

Only for the first 4 of these 8 Q- values does a P exist such that D X P = Q ■ (u" + u), 
i.e. only 4 first integrals are found, the remaining 4 solutions represent only adjoined 
symmetries. This inability to find 5 first integrals is not a weakness of the computer 
program but the price to pay for the fourth method to have shorter conditions (fT^j) 
compared with the conditions (|13|) of the third method. If the conditions (J 14)) of the 
fourth method are less restrictive than the conditions (|13|) of the third method, why then 
does it find fewer first integrals, 4 instead of 5? The answer is that the 5 th conservation 
law is contained in the second half of the 8 solutions: the 8 th solution for Q plus 2 times 
the 7 th solution gives a value for Q that is not only an adjoined symmetry but also an 
integrating factor for an additional first integral. 
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To summarise, the fourth method gives shorter, more manageable conditions which, 
strictly speaking, do not compute conservation laws but adjoined symmetries. Conser- 
vation laws can be derived through appropriate linear combinations of adjoined sym- 
metries which has to be investigated separately. This theoretical weakness of the fourth 
method does usually play no role in practical applications. 

Applying the program ConLaw2 that corresponds to the above method to the sine-Gordon 
equation (fTTj) gives the following table. The typical feature of this approach is the slower 
increase of the size of conditions. Apart from the order case they increase by a factor 
of about 2 which itself is increasing slightly with the order. Compared with the previous 
method the size of conditions grows slower which allows to go higher in the order. Because 
the conditions that are generated are only necessary, not sufficient, they are slightly more 
difficult and expensive to solve. This causes longer running times for low order investigations. 
Time limitations could be overcome to some extend by faster computers. 
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order 
of Q 


r 

no oi 
terms 


independent 
variables, [no of var.] 


functions to 
compute, [no of arg.] 


cons, 
laws 
found 


time to 
solve dHJ) 
h:min:sec 





6 


[5] 


Q(t,x,u), [3] 


none 


1 sec 


1 


21 


t,x,u {2) /utx, [7] 


Q{t, xM% [5] 




4.3 sec 


2 


45 


t,x,u®/u£\ [9] 


Q(t,x,n( 2 )/n te ), [7] 


none 


12 sec 


3 


99 


t,x,u^/u t i 2 \ [11] 


Q(t,x,u®/u£ } ), [9] 


GD 


50 sec 


4 


202 


t,x,u®/u£\ [13] 


Q{t,x,u^/u^), [11] 


none 


2 min 43 sec 


5 


435 


t,x,u^/ Ut { *\ [15] 


Q(t,a;,n( 5 )/ni 3) ), [13] 




16 min 10 sec 


6 


870 


t,x,u^/u t f \ [17] 


Q(t,x,u^/u t ^), [15] 


none 


49 min 20 sec 


7 


1836 


t,x, M ( 8 )/ M i 6) , [19] 


Q{t,x,u^ /u^), [17] 


dsnj) 


8 h 


8 


3643 


f,z, [21] 


Q(t,x,vF>/u£>), [19] 


none 


5 h 22 min 


9 


7434 


t,x,uW/u£\ [23] 


Q(t,a;,n( 9 )/ni 7) ), [21] 


dSH) 


25 h 



Table 4: The program ConLaw2 applied to compute conservation laws of the 
sine-Gordon equation. 
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4.5 Overview 



The circumstance that the number of conservation laws for the sine-Gordon equation that 
were found by the different methods varies is due to the varying computational complexity of 
the determining equations they generate. The first three methods will find all conservation 
laws if memory and time requirements would not matter. The fourth approach is different in 
that it generates only necessary conditions which are often sufficient (if they have the same 
solution set as the other methods) but sometimes not. In that case the conditions (jHj) have 
additional adjoined symmetries as solutions. It may be that only specific linear combinations 
of them give a conservation law, as demonstrated with the example in the previous section. 
The following comments concentrate on complexity issues and other characteristic differences 
between the four approaches. 

Arranging the methods as in the table below one can compare rows I, II and columns A,B. 





A 


B 


I 


DivP| A=0 = 


Xv,A-D)j (qA) =0 Vi/ 

\ J / A=0 


II 




£j(-0)j(^£ m Q,A.)=O Vz/ 



Table 5: The four approaches arranged in a table. 



I-II: The conditions in row I are to be solved in the space of solutions (|a=o), i n row II they 
are not. This means that methods of row I can not be applied if equations or constraints 
A M = can not be solved for a leading derivative but methods of row II can. Due to these 
substitutions the conditions in row I have fewer terms and involve fewer different derivatives 
of u than conditions in row II. The complexity of conditions and the number of conservation 
laws up to some order obtained in row I depend on whether A^ = is used to substitute 
lower order u- derivatives by higher ones or higher ones by lower ones. There are two reasons 
for this. 

1) Substitutions based on = A in Q may give extra restrictions for Q. For example, 
determining Q for conservation laws of the Korteweg de Vries equation = A 
and restricting Q to be of 2 nd order, then a substitution u t = u xxx + uu x would imply 
Q = Q(t,x,u,u x ,u xx ), whereas a substitution u xxx = u t — uu x would not restrict Q = 
Q(t, x, u, u t , u x , u tt , u tX} u xx ). 
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2) If a lower n-derivative is substituted by higher ones using = A in the conservation law 
conditions in row I then such substitutions may increase the order of w-derivatives in which 
the conservation law conditions have to be satisfied identically. By that the desired effect of 
lowering the number of ^-derivatives in which the conditions have to be fulfilled identically 
is lost. For example, condition IB for A = u a — u xxt is = Q tt + ^{Qu xxt ) xx t which includes 
up to 6 th order w-derivatives (if Q is not of higher than 3 rd order). By substituting u tt = u xxt 
the order would increase to seven. 

Hence, substituting lower order ^-derivatives by higher order w- derivatives gives more over 
determined conditions for a less general ansatz. Such conditions are easier to solve, which 
may allow higher orders of Q to be investigated. However, one then may miss conservation 
laws of some order in P or Q. 

These aspects are not an issue in row II as no substitutions are made there. 
A-B: In column A the single first order conservation law condition itself is to be solved, 
and in column B the integrability conditions of column A, which result when the conserved 
current P is eliminated are to be solved. Conditions in column B involve as many equations 
as there are functions and they are of the same order as the highest derivatives of in 
A„ = 0. 

The methods in column A compute P and therefore also trivial conservation laws when 
determining the general solution of the determining equations. This may complicate the 
solution of the determining equations to some extend. After the determining equations are 
solved the trivial conservation laws are easily dropped. For method IA one simply checks 
whether Div P = holds identically and in method IIA trivial conservation laws have zero 
integrating factors Q^. The general solution of the determining conditions in column B do not 
generate trivial conservation laws. Also, conditions in column B are more straightforward 
to solve, they can be separated with respect to many high order u-derivatives and yield 
highly over determined systems. The disadvantage of methods in column B is that already 
their formulation may exceed available computational resources. Another potential problem 
with using methods in column B is the following. If one or more linear PDEs from the 
over determined conditions remain unsolved (for example, when investigating the Burgers 
equation (|T7|) then the heat equation remained unsolved (see equations (|T7|) .(fTH j) ) then the 
program will usually not be able to compute P l from the Q^. A way out is to use methods 
IB or IIB to get and to use that as input to get P l from method IIA. 

Differences between the approaches are amplified with problems that involve an increasing 
number of PDEs and an increasing number of independent variables. A recommendation for 
tackling a new single PDE / system of PDEs would be: 

• If any ansatz for some or all of the P l is to be made then the method IA has to be used. 
(An example is the question whether an easily integrable conservation law D X P X = 
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with P y = exists.) 

• Otherwise try IIB first. 

• If the conditions become too large to handle then try IB or even IIA or IA. 

• If IIB or IB provide the but not the P % (for example because arbitrary functions 
appear or some conditions remain unsolved) then try IIA using the computed as 
input. 

5 The computer algebra programs 

The names of computer algebra programs for the four approaches are: IA: ConLaw1,IB: 
ConLaw2, IIA: ConLaw3, IIB: ConLaw4. They and the program Crack for solving the 
over determined conditions are written in the computer algebra system Reduce (ConLawI 
through ConLaw4 by the author, Crack by the author and Andreas Brand). From the 
general solution individual conservation laws are extracted by picking one arbitrary constant 
or function and setting all other arbitrary constants or functions to zero. The problem is to 
find, whether all arbitrary constants and functions are independent or whether some can be 
dropped without loss of generality. This problem itself leads to an over determined system of 
conditions which in general is very over determined and easy to solve. A description of that 
method is given in [SS] where also the computation of Q from P and P from Q is explained. 

Compared to other computer algebra programs, the package Crack has a wide variety of 
techniques for solving over determined PDE-systems. This allows the following new features 
as compared with other computer programs, a list of which and a short description is given 
in H2|: 

• In all four computer programs P as well as Q is computed. 

• By solving systems of over determined differential equations it is possible to find con- 
servation laws with non-polynomial, even non-rational P, Q. 

• If memory requirements are not too high then the program will make a definite state- 
ment about the existence of conservation laws of a given order. In the majority of these 
cases the program will find the explicit form of the conservation law, otherwise it will 
return unsolved equations. 

• It is possible to find conservation laws with an explicit dependence of P, Q on the 
independent variables. 
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• There is no limit on the number of DEs nor the number of independent variables to 
be investigated for conservation laws other than a limit through the complexity of 
computations. Although not demonstrated in this paper, the program is able to handle 
ordinary differential equations (ODEs) as well. 

• It is possible to determine values of parameters in the DE such that conservation laws 
exist. 

• For each of the four programs ConLawI through ConLaw4 an ansatz for P l and/or 

can be input to specify to some extend conservation laws to be calculated. 

A program written by G6kta§ and Hereman jS] for computing conservation laws of PDEs 
makes a polynomial ansatz for conservation laws and finds the coefficients in this ansatz by 
solving a linear algebraic system of equations. Compared with that, the programs ConLawI 
through ConLaw4 are able to find more general conservation laws and to make a definitive 
statement in case the order is not too high to complete the computations. On the other 
hand, the program of G6kta§ and Hereman was later extended to handle differential- difference 
systems , HD| , HH . 

Before showing examples which highlight the special abilities of ConLawI through Con- 
Law4 a comment to the treatment of ODEs will be made. Although all methods and pro- 
grams are applicable equally well to ODEs, the form of the ansatz for the integrating factor 
Q or for the first integral P to be made will usually be different. An n th -order ODE has 
always first integrals of order n — 1 and any arbitrary function of first integrals is a first 
integral as well. In order to obtain an over determined system of conditions, the ansatz for 
a first integral must not contain functions of all variables x,u,u', . . . ,d n ~ 1 u/dx n ~ 1 but, for 
example, a polynomial in d n ~ 1 u/ ' dx n ~ x with arbitrary functions of x,u,u', . . . ,d n ~ 2 u/dx n ~ 2 
as coefficients or any other combination of functions of less than n+1 variables, see also [I] 
for more details. Special features of the ConLaw programs that are not available with other 
programs are demonstrated with the following examples. 

Example: The Ito equations for two functions u = u(t,x), v = v(t,x) read 

Ut = u xxx + Quu x + 2vv x 
v t = 2(uv) x . 

The first 7 conservation laws calculated by the program ConLawI which in turn calls Crack 
to solve condition (|1U|) . have the following values of P l : 

u, v, u 2 + v 2 , u 2 x - 2u 3 - 2uv 2 , (Auv 2 - v 2 x ) / v 3 , 

u 2 x - lOuul - ^v x u x + 5n 4 + 6u 2 v 2 + v 4 , {{2vv xx - Auv 2 - 3v 2 ) 2 + 16v 6 ) /v 7 . 
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(P x is not shown due to its length. It could be computed easily from P f .) What is interesting 
is that 2 of the 7 conservation laws have a non-polynomial expression for P l and as far as 
the author knows these conservation laws have not been known so far. 

Example: The following equations [THj describe low-frequency Alfven waves propagating 
parallel to an external magnetic field in a relativistic electron-positron plasma [HI]. Typical 
for them is the symmetry with respect to interchanging the two functions u = u(t,x),v = 
v(t, x) due to the same charge-to-mass ratio for both kinds of particles. The equations are 



A 2 



u t + r x 
v t + s x 



0. 
0. 



with 
with 



r = u(u 2 + v 2 ) + u x 
s = v{u 2 + v 2 ) + v x 



(15) 



The equations themselves have the form of conservation laws. We find the following additional 
ones: 



AuA l +AvA 2 = D t [2(u 2 + v 2 )} + 

D x [Auu xx - 2u 2 + Avv xx - 2v 2 + 3(u 2 + v 2 ) 2 } 
4rAi+4sA 2 = D t [(u 2 + v 2 ) 2 - 2u 2 - 2v 2 } + 

D x [Au t u x + Av t v x + 2u xx + 2v xx + A(u 2 + v 2 ) x 
((3(u 2 + v 2 )t - x)(uu t + vv t ) + uu xx + vv xx "j 
A(xu - 3tr)Ai + A(xv - 3ts)A 2 = 



D, 



D x 2 



+ 



3t((2u 2 + 2v 2 - (u 2 + v 2 ) 2 ) + 2x(u 2 + v 
(uu t + vv t ) (-x 2 + (u 2 + v 2 ) (6tx - 9t 2 (u 2 + v 
3t(u 2 + v 2 f + 3x(u 2 + v 2 ) 2 + 2x{uu xx + vv xx ) - 3tr 
■3ts 2 — 2uu x — 2vv x — xu 2 — xv 2 — 6tu t u x — 6tv t v x 



Whereas the first two are known [SI] , the last one shows an explicit x, t-dependence and is 
new. Further investigation provides that no conservation laws exist with the characteristic 
functions of 3 rd or 4 th order (if u t ,v t are substituted using (JT3|) ) . 
Example: The following equation of Gibbons and Tsarev [7] 







u, 



+ UyU 



y Uj %y 



(16) 



for u = u(x, y) is unusual in that it has already 5 conservation laws of first order. The 
characteristic functions contain x, y explicitly. Up to first order they are: 



1. 



3u + 2u x + 3x, 2u + 3u x u y + Au y x + y 



22 



10u y + 6u 2 x + 2Au x u 2 y + 20u x x + 30u 2 y x + I2u y y + 2u + I5x 2 , 
3w° + 6u x u y + 10u x u y + 18u x u y x + 4:U x y + 12u y x + 6w y ?/ + 12u y x + 2u y u + 6xy. 
Example: The Liouville equation for a function w = u(x,y) reads 

A = u xy - e u . 

Conservation laws of order zero found by ConLaw2 are 

(f x + fu x )A = D x (-e u f)+D y (f x u x + fu 2 /2), f = f(x) arbitrary 

(g y + gu y )A = D y (-e u g) + D x (g y u y + gu 2 /2), g = g(y) arbitrary. 

Because the ansatz made is investigated in full generality, any free functions in the conserva- 
tion law will be found if the conditions can be solved completely by Crack. Otherwise the 
remaining conditions are returned as in the following example. 
Example: The Burgers equation in the form 

A = u t - u xx - -ul = 0, u = u(t,x) (17) 

has zeroth order conservation laws 

fe u ' 2 A = D t (2fe u ' 2 ) + D x (e u ' 2 (2f x - fu x )) (18) 

with f = f(t, x) satisfying the linear reverse heat equation = f t + f xx . 8 

The occurrence of free functions in the conservation law indicates linearizability of A = 0, 
which is the case for both previous examples. The following example involves more than 2 
variables. 

Example: The Kadomtsev-Petviashvili equation for u = u(t,x,y) with the abbreviation 

w = Ut + 2uu x + u xxx 

is 

= A = W X - Uyy. (19) 

8 Although already used in j^j, [331 this example is shown again as it also serves to demonstrate an extension 
to non-local conservation laws in section |5| 
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Its zeroth order conservation laws include an arbitrary function c = c(t): 



(6cxy + Q|/ 3 )A 



cyA 

{2cx + c t y 2 )A 



cA 



D X (CW) + Dy(-CUy) 

D x (cyw) + D y (cu - cyu y ) 
A(-2cu) + 

D x \i2cx + c t y 2 )w — leu xx — 2cu 2 ) + 
L> y (-(2ct + c t y 2 )uy + 2qw?/) 
D t (-Qcyu) + 

((6cx?/ + Q?/ 3 )w - eq/Ua-e - 6q/w 2 ) + 

-D y ( — (6cxy + c t y 3 )u y + 3c t uy 2 + Qcxu) . 



(20) 
(21) 



(22) 



(23) 



It is somewhat remarkable that although equation (|19j) does not involve Ut but only u x t 
nevertheless the conserved density P t in the last two conservation laws involves u and not 



In the following section we give examples for an extension of our method to compute non- 
local conservation laws and report on the possibility to determine parameters in the equation 
such that conservation laws exist. 

6 Extending applicability 
6.1 Non-local conservation laws 

The implementations of the four methods have a common limitation: the characteristic 
functions Q and the conserved current P must depend functionally only on a finite number 
of derivatives of the u a . No dependencies on integrals are possible. The same restriction is 
usually made when generators of Lie-symmetries are determined for differential equations. 
Whereas this restriction is less severe when calculating symmetries of PDEs, it is a serious 
restriction for the determination of conservation laws. To give an example, Burgers' equation 
in the form 



has as low order conservation law only the trivial one D t u — D x (u x + u 2 /2) = 0. In order to 
include dependencies on / udx one could set u = v x for some function v(x,t) and investigate 
conservation laws depending on v and derivatives of v. For Burgers' equation such a substi- 
tution alone is not enough. In addition one has to realize that (|24|) can be integrated with 





(24) 
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respect to x to f(t) t — v t — v xx — v 2 for some function / = /(£). Renaming v — / — > u gives 
(fT7|) and its conservation laws (JTSJ). 

To give a further example, we consider the Boussinesq equation describing surface water 
waves whose horizontal scale is much larger than the depth of the water [T].|13| 

Utt - u xx + 3uu xx + 3u x + au xxxx = 0. (25) 

Calculating conservation laws, using (|23j) to substitute u xxxx , the only characteristic functions 
Q up to 4 th order are l,x,t, xt. On the other hand, substituting u = v x , integrating (|2*3j) with 
respect to x and renaming v — f — > v gives 

vtt - Vxx + 3v x v xx + av xxxx = (26) 

having 2 conservation laws with characteristic functions l,t which x-differentiated give the 
conservation laws above with characteristic functions 1, t. In addition two new conserva- 
tion laws with characteristic functions v x ,v t result. Repeating this step again: v = w x , 
x-integration of (J2HJ), w — f — >• w gives 

w a -w xx + 3/2w xx + aw xxxx = (27) 

with three third-order conservation laws. Two of them have characteristic functions w xxx , w xx t 
which correspond to the above conservation laws with characteristic functions v x ,v t . In 
addition one extra conservation law with Q = w txx — w txx w xx + w tx w xxx — ^w t tt exists. 

A third example is the Kadomtsev-Petviashvili equation already discussed above. 9 After 
a substitution u — v x , x-integration of (|T^|) and v — f — > v the equation is 

= [v t + v xxx + vl\ x - Vyy. 

Apart from conservation laws equivalent to ([20 )1 .([21 )1 three new conservation laws result with 
characteristic functions 

-c tt y 2 - 2c t x + Acv x 

-c 3 ty 3 - Gcuxy + I2c t yv x + 24cv y , 

-c it y 4 - 12c 3t xy 2 + 2Ac tt y 2 v x - I2c tt x 2 + 48c t xv x + 96c t yv y + A8c t v + \AAcv t . 

Conserved currents are omitted due to their length. Repeating this transformation again 
does not yield conservation laws with characteristic functions of order less than three. 

The purpose of this paragraph was to show that even if computer algebra programs 
ConLaw, Crack do only allow the investigation of local conservation laws depending on 
a finite number of derivatives of the unknown functions, we still may be able to enlarge the 
range of search by a contact transformation and integration of the PDE. 

In the next section we extend the computation of conservation laws to the computation 
of parameters such that conservation laws exist. 

9 The hint to try KP for this extension was given by Alan Fordy. 
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6.2 Differential equations with parameters 

In applications it is common that the DEs contain parameters and usually it would be 
desirable to know conservation laws which are valid for all possible values of these parameters. 
But as the example below shows, often conservation laws exist only for special values of 
parameters. Even if these parameter values are not of interest from the application side of 
view, the conservation laws valid for these values can at least be used, for example, to test 
numerical code. Another purpose for determining parameters together with conservation 
laws could be to find integrable equations from a more general class of equations. 

The problem to determine parameters such that conservation laws exist is potentially 
much harder than determining conservation laws which are valid for any values of these 
parameters. This is because the conservation law determining equations become non-linear. 
Expressions may become unmanageably large and many sub cases may have to be considered. 
To use ConLawI through ConLaw4 for such calculations one only has to specify in its call 
the names of parameters to be computed (more details in the ConLaw manual). 

Example: The 5 th -order Korteweg - de Vries equation 

u t + au 2 u x + (3u x u xx + juu 3x + u 5x = (28) 

with constant parameters a,/3,7 includes well known special cases 0,0, |Hj, [IE], [21] : for 
a = 30, (3 = 20, 7 = 10 the Lax equation [TH], for a = 5, (3 = 5, 7 = 5 an equation due to 
Sawata, Kotera [2H] and Dodd and Gibbon [5], for a = 20, (3 = 25, 7 = 10 an equation due 
to Kaup [T7j and Kupershmidt, for a = 2, (3 = 6, 7 = 3 an equation due to Ito [To] . 

The following zeroth and first order conservation laws are calculated with ConLawI 
(omitting P x due to its length in the last two of these conservation laws): 

• Q = l, P* = u 

• a = /5 = 7 = 0: As (|2*Hj) becomes linear, a conservation law is obtained with a 
characteristic function Q = Q(x,t) satisfying the adjoint PDE Q t + Q 5x = with 
P* = Qu. 

m a = 0, 7 = P/3 : Q = x 2 , P l = x 2 u 

• a = 0, 7 = (3/3 : Q = x, P t = xu 

• 7 = ^/2: Q = 2u, P l = u 2 

• a = ^(-2(3 2 + 7/5 7 - 3 7 2 ): 

Q = Q0u xx t((3 - 37) + Qu 2 t(2(3 2 - 7/?7 + 37 2 ) + QOx 
pi = 30u 2 x t{-(3 + 37) + uH{Af3 2 - 14^7 + 67 2 ) + 6Qux 
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. « = ^(-2/3 2 + 7/3 7 - 3 7 2 ): 
Q = 30u xx + 3u 2 (2f3 -7), 
pt = _15 U 2 +U 3( 2( 5_ T ) 

We find the same conservation laws as found by the program of Gokta§ and Hereman and in 
addition a few conservation laws with explicit x, t-dependence. 

7 Summary 

Four approaches to find conservation laws have been compared with respect to their com- 
plexity and other characteristic features. 

In a number of examples, conservation laws have been given, some of them new, which 
show that the programs ConLaw1,...,ConLaw4 and Crack can be used to find local, 
not necessarily polynomial, conservation laws with explicit variable dependence and free 
functions. The programs are, in principle, applicable to problems with arbitrarily many 
equations, functions and variables. 

The package is part of the Reduce distribution. The latest version of the programs 
including a manual and a test file are available from 
http : //lie .math. brocku. ca/twolf /crack/ crack. tar . gz . 
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9 Appendix: 

Conservation Laws of the sine-Gordon equation 

In this appendix conservation laws for the sine-Gordon equation 

u tx - sin(u) = 

are shown as they have been computed by ConLawI-4 and as they are refered to in tables 
above. They are not new, we provide them only to illustrate computer results. Except for 
the first conservation law, for all the following there is an additional conservation law due to 
the x <-> t symmetry. These results are further examples of the ability of the programs to 
compute non-polynomial conservation laws. 

2{tu t — xu x )(u tx — sin(w)) = D t 2cos(u)t — u x x + D x —2cos(u)x + uft (29) 
2u t (u tx - sin(u)) = D t [2 cos(u)] + D x [uf\ (30) 
(8u 3t + 4ut)(u tx - sin(u)) = D t 4 cos(u)u^ + 8u tx u u - Su tt sin(u) + D x -Au 2 tt + u\ (31) 
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2(-8u 5t - 20u 3t u 2 - 20u 2 t u t - 3uf)(u tx - sm(u)) 
D t 2 —8 cos(u)u st u t + 4 cos(u)u 2 t — 3 cos(u)u^ — 8u tx u 4t — 20u tx u tt u 2 
+8uu sin(w) + 8u 3t u ttx + I2u tt u 2 sin(u) 



(32) 



+ D x -8u 2 3t + 20u 2 tt u 



8(—lQu 7t — 56u 5t u t — 22AuuUttUt — lQ8u 3t u t — 280u 3t u tt — 70u 3t u t — 140u tt u t 
x(u tx - sm(u)) 



5ui 



(33) 



A8 



16u 6t sm(u) — 16u tx u 6t — 16 cos(u)u 5t u t + 16 cos(u)u4 t u tt — 8 cos(u)u 



— 40 cos(u)u 3t u t — 20 cos(u)u tt u t — 5cos(u)u t — 56u tx U4 t u t — I12u tx u 3 tu tt ut 
-56u tx u% - 7§u tx u tt u\ + 16u 5t u t tx - lQu 4 tU 3tx + 40u 4t u 2 sin(u) + 56u 3t u ttx ul 
+ 160u 3t u tt u t sin(u) + 40w^ sin(-u) + 30u tt u 4 sin(w) 



+ D x 



64m 



At 



224u 2 3t u 2 + \\2u% + 28§u\ t u\ 



5m 



2(— 128-ugt — 576u 7t u 2 — 3456u 6t u tt u t — 7296u 5t u 3t u t — 6720u 5t u 2 t — 1008«5 t u^ 



-44\<ou\ t u t - 24l92u 4t u 3t u tt - 8064u 4t u tt u 3 t - 5824i4 
-840u 3t u 6 t - 6384u 4 tt u t - 2520u 2 t u 5 t - 35uf) 
x(u tx - sin(u)) 



6048^X 



24864u 3t u 2 t u 2 



= D t 2 



(34) 

128-Ugt sin(u) — 128 cos(u)u 7t u t + 128 cos(u)uQ t u tt — 128 cos(u)u 5t u 3t + 64 cos^)^ 

— 448 cos(u)u5 t Uf — 1344 cos(u)u 4t UuU 2 — 1568 cos{u)u\ t u 2 — 1344 cos(u)u 3t u 2 t u t 
—560 cos{u)u 3t u\ + 336 cos{u)u tt — 840 cos(u)u 2 t u 4 — 35 cos(w)w® — 1008u tjX U4 t u t 
-576u t , x u 6t u 2 - 230Au tjX u 5t u tt Ut - 4992u t , x u At u 3 tUt - 4416u tiX u At u 2 t - 128u tjX u 8 t 
-518Au tjX ul t u tt - 4032u t:X u 3t u tt u 3 - 4256u tjX u%u 2 - 840u t , x u tt Uf + I28u 7t u tt , x 
-128u et u 3tjX + 448u 6t u 2 sm(u) + 128u 5t u 4t>x + 576u 5t u tt<x u 2 + 2688u 5t u tt u t sm(u) 
-576u At u 3t)X u 2 + 4480u4tU 3 tUt sin(u) + U52u4tu tt , x u tt u t + 4032u 4t u 2 t sin(u) 
+560-U4 t -u t sin(-u) + 1920u 3t u t t, x u t + 5824-u 3t -u tt sin(-u) + 3264u 3t Uu, x u tt 
+1008u 3t u t t,xUt + 4480u 3 tU U Ut sm(u) + 3360u^ t u 2 sm(u) + 280u tt Ut sm(u) 
+ D x [-128u 2 5t + 576u 2 4t u 2 - 1280ul t u t - 3264u 2 3t u 2 tt - 1008x4*4 + 2128u Wt 
+840u 2 tt u 6 t - 7ul° 
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